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Abstract In a series of papers we have recently studied the clustering of LRG 
galaxies in the latest spectroscopic SDSS data release, which has 75000 LRG galax- 
ies sampling 1.1 Gpc^^/h^ to z=0.47. Here we focus on detecting a local maxima 
shaped as a circular ring in the bidimensional galaxy correlation function ^ {a, n), 
separated in perpendicular a and line-of-sight n distances. We find a significant de- 
tection of such a peak at r ~ 1 lOMpc/h. The overall shape and location of the ring 
is consistent with it originating from the recombination-epoch baryon acoustic os- 
cillations (BAO). This agreement provides support for the current understanding of 
how large scale structure forms in the universe. We study the significance of such 
feature using large mock galaxy simulations to provide accurate errorbars. 



1 Gravitational instability 

Is the large scale structure that we see in the galaxy distribution produced by grav- 
itational growth from some small initial fluctuations? We will explore two ways of 
addressing this question with measurements of the 2-point galaxy correlation: 

^(r,r) = (5(n,05(r2,r)) (1) 

where r = | r2 — r 1 1 and 5 (r) = p (r) /p — 1 is the local density fluctuation about the 
mean p = (p), and the expectation values are taken over different realizations of the 
model or physical process. In practice, the expectation value is over different spatial 
regions in our Universe, which are assumed to be a fair sample of possible realiza- 
tions (see Peebles 1980). The measured redshift distance of a galaxy differs from 
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the true radial distance by its peculiar velocity along the line-of-sight. We can split 
the distance r into its component along the line-of-sight (LOS) % and perpendicular 
to the LOS C7, where = + o^. Azimuthal symmetry implies ^ is in general a 
function of K and a alone: ^ (cj, n). 

Consider the fully non-linear fluid equations that determine the gravitational evo- 
lution of density fluctuations, 5, and the divergence of the velocity field, 0, in an ex- 
panding universe for a pressureless irrotational fluid. In Fourier space (see Eq. 37-38 
inffl): 

5 + e = - Jdkidk2 a{hM)Q{h)d{k2) (2) 
d + J^d + lQ.,„H-8 = - Jdkidk2 P{ki,k2)0{ki)e{k2) 

where derivatives are over conformal time dt = dt/a and Jf{T) = dlna/dT = aH 
is given by the expansion rate H = d/a of the cosmological scale factor a. On the 
left hand side 5 = 5{k, t) and 9 = 9{k, t) are functions of the Fourier wave vector 
k. The integrals are over vectors k\ and k2 constrained to ^ = ^9 ^ ■ The right hand 
side of the equation include the non-linear terms which are quadratic in the field and 
contain the mode coupling functions a and j3 . When fluctuations ai-e small we can 
neglect the quadratic terms in the equations and we then obtain the linear solution 
5l- The first equation yields 8l = —9l, which combined with the second equation 
yields the well known harmonic oscillator equation for the linear growth: 

dL + .yfdL-^n,n,yf^dL^O (3) 

Because the Fourier transformation is linear, this equation is valid in Fourier or in 
configuration space. In linear theory each Fourier mode and each local fluctuation 
evolves independently of the others, moreover, they all grow linearly out of the 
initial fields with the same growth function, ie ^/.(f ) = D{t)5o, where 80 is the value 
of the field at a point (or a given Fourier mode) at some initial time and D{t) is 
the linear growth function, which is a solution to the above harmonic equation. In 
a flat universe dominated by cold dark matter (CDM) we have that H <^ ot'^I'^ and 
Z)(f ) goes as the scale factor D(t) oc a. In an accelerated phase, such as ACDM, the 
growth halts or grows less rapidly with a. Thus measurements of D{t) can be used 
as an independent diagnostics for accelerated expansion. 



1.1 BAO signature 

Consider now our observable, the 2-point function in Eq.ITl. In linear theory: 
^ (r,f ) =D(f)^^ ('"jO)- This means that on large scales, ie r > lOMpc/li, where 5 < 1, 
we then expect the shape of £, (r) today to be the same as the shape it took in the 
early universe. The prediction, ignoring redshift space distortions, is shown in Fig[T] 
The BAO ring corresponds to the local maxima at a radius of about 1 IQMpc/h. This 
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Fig. 1 The prediction of 
^{a,7t) without redshift space 
distortions. The vertical axis 
shows the radial direction, 
7t, while the horizontal panel 
shows the transverse direction 
CT. There is a prominent local 
maximum corresponding to 
the BAO ring at a radius 
of about llOMpc/h and 
~ IQMpc/h thickness (wide 
green circle between blue 
rings). 
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will be tested here by comparing the linear theory prediction for the baryon acoustic 
oscillation (BAO) with the measurements in the large scale galaxy distribution. 

The mean BAO signature in the 2-point correlation has been detected in Lumi- 
nous Red Galaxies (LRG's) of the SDSS galaxy survey Q using the monopole, ie 
the average signal of ^ (tt, ff) in circles of constant TT^ + ff^. In |l3]|4l and below we 
will show that the galaxy distribution has a BAO ring very similar to that predicted 
in the initial conditions of the CDM model. 



1.2 Growth factor 

On the other hand, if we knew the shape and amplitude of the initial conditions 
B, (r, 0), we could then estimate D(f )^ ~ ^ (r, f )/^ (r, 0) from measurements of ^ (r, f ) 
and compare it to the linear solution of Eq.fS) to test gravitational instability. But 
this approach is difficult in practice because there is a bias in the amplitude of 
galaxy clustering compare to the one in dark matter fluctuations. We could in- 
stead test gravitational growth with independence of time or initial conditions by 
using the linear relation between density and velocities, 5i = —Oi = fSi, where 
/ = D/D = dlnD/dlna is the velocity growth factor For a flat dark matter domi- 
nated universe D^a and / = 1 , while for a flat accelerated universe / = Q,„ (z)"^ < 1 , 
where 7 is the gravitational growth index and I2„,(z) is the matter density at a red- 
shift z where a — 1/(1 + The growth index separates out two physical effects on 
the growth of structure: Q,„ (z) depends on the expansion history while 7 depends 
on the underlaying theory of gravity Q. The value 7 = 0.55 corresponds to stan- 
dard gravity, while 7 is different for modified gravity, for example 7 = 0.68 in the 
braneworld cosmology. We will show next how / can be measured using redshift 
space distortions in the galaxy correlation function ^ (r). 
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In this work we use the most recent spectroscopic SDSS data release, DR6 (Q). 
We use the same samples and methodology here as presented in fS) of this series. 
LRG's are targeted in the photometric catalog, via cuts in the (g-r, r-i, r) color-color- 
magnitude cube. We need to k-correct the magnitudes in order to obtain the absolute 
magnitudes and eliminate the brightest and dimmest galaxies. We have seen that the 
previous cuts limit the intrinsic luminosity to a range —23.2 <Mr< —21 .2, and we 
only eliminate from the catalog some few galaxies that lay out of the limits. Once 
we have eliminated these extreme galaxies, we still do not have a volume limited 
sample at high redshift. For the 2-point function analysis we account for this using 
a random catalog with identical selection function but 20 times denser (to avoid 
shot-noise) . The same is done in simulations. 

There are about 75,000 LRG galaxies with spectroscopic redshifts in the range 
z = 0.15 — 0.47 over 13% of the sky. We break the full sample into 3 independent 
subsamples with similar number of galaxies; low z = 0. 15 — 0.30, middle z = 0.30 — 
0.40 and high z = 0.40 — 0.47. In this analysis we will just show results for the 
z = 0. 15 — 0.30 sample (see |3]|4J for other samples). 

To estimate the correlation £, {a, 7z), we use the estimator of [Sl, 

, DD 2DR+RR 

^(^.^)- m 

with a random catalog Nr = 20 times denser than the SDSS catalog. The random 
catalog has the same redshift (radial) distribution as the data, but smoothed with 
a bin dz = 0.01 to avoid the elimination of intrinsic correlations in the data. The 
random catalog also has the same mask. We count the pairs in bins of separation 
along the line-of-sight (LOS), n, and across the sky, a. The LOS distance n is just 
the difference between the radial comoving distances in the pair. The transverse 
distance a is given by x/i^ — tt^, where s is the net distance between the pair. We 
use the small-angle approximation, as if we had the catalog at an infinite distance, 
which is accurate until the angle that separates the galaxy pair in the sky is larger 
than about 10 degrees for £, {a, n) (see IH and ifTOl ). This condition corresponds to 
transverse scales larger than a = %QMpc/h for our mean catalog. 

The right panel in Fig|2] shows the measurements of i^(;r,CJ) in the z = 0.15 — 
0.30 sample. As we will show below there is a remarkable agreement with the pre- 
dictions (left panel) and there is good evidence for a BAO ring. 



2.1 Errors and simulations 

There are two sources of error or variance in the estimation of the two-point corre- 
lation: a) shot-noise which scales as one over the square root of the number of pairs 
in each separation bin b) sampling variance which scales with the amplitude of the 
correlation. It is easy to check that for the size and number density of our sample. 
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Fig. 2 A comparison of ^{a.Tl) in data and models for the z=0. 15-0.30 slice. The vertical axis 
shows the radial direction, 7t, while the horizontal panel shows the transverse direction a. The left 
panel includes a model of redshift space distortions which gives the best fit to the monopole in 
the galaxy data (see \3 \). Right panel shows the LRG SDSS measurements using the same color 
scheme. The data shows a prominent BAO ring at a radius of about 1 iOMpc/h, in good agreement 
with the model. 



the shot-noise term dominates over the sampUng variance error This has been con- 
firmed in detail by using numerical simulations (see 13] for more details). The simu- 
lation contains 2048^ dark matter particles, in a cube of side 76S0Mpc/h (which we 
call MICE7680), Qm = 0.25, Qh = 0.044, CTg = 0.8, = 0.95 and h = 0.7. We have 
divided this big cube in 3^ cubes of side 2 x 1275Mpc/h, and taking the center of 
these secondary cubes as the observation point (as if we were at z=0), we apply the 
selection function of LRG, which arrives to z=0.47 (r=1275Mpc/h). We can obtain 
8 octants from the secondary sphere included in the cube, so at the end we have 8 
mock LRG catalogs from each secondary cube, which have the same density per 
pixel as LRG in order to have the same level of shot noise, and the area is slightly 
smaller (LRG occupies 1/7 of the sky with a different shape). The final number of 
independent mock catalogs is M = 216 (27x8). We also apply redshift distortions in 
the line-of-sight direction.? = r + Vr/H{z)/a{z), using the peculiar velocities Vr from 
the simulations. The error covariance is found from the dispersion of M realizations: 

1 M 

cu = TfLi^iif-mmjf-^{j)) (5) 
™ k=i 

where ^{i)'' is the measure in the k-th simulation (k=l,...M) and ^{i) is the mean 
over M realizations (which we have checked that agrees with overall mean, indicat- 
ing that volume effects are small). The case i=j gives the diagonal error (variance). In 
our analysis we model the errors using dark matter groups. These groups are chosen 
to have the same number density and amplitude of clustering as the observed LRG's. 
The resulting errorbars from simulations are typically in good agreement with Jack- 
knife errors from the actual data (see lO for details). We have also checked in fS) 
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that we can recover the theory predictions for ^{a,n) within the errors by using 
mock simulations with similar size as the real data. 



2.2 Redshift space distortions 

Radial displacements caused by peculiar velocities lead to redshift distortions, with 
two important contributions. The first, on large scale fluctuations, caused by coher- 
ent bulk motion. We see walls denser and voids bigger and emptier, with a squashing 
effect in the 2-point correlation function along the line-of-sight: known as the Kaiser 
||6l effect. At small scales, random velocities inside clusters and groups of galaxies 
produce a radial stretching pointed at the observer, known as fingers of God (FOG). 
Although such distortions complicate the interpretation of redshift maps as posi- 
tional maps, they have the advantage of bearing unique information about the dy- 
namics of galaxies. In particular, the amplitude of distortions on large scales yields 
a measure of the linear redshift distortion parameter /. 

In the large-scale linear regime, and in the plane-parallel approximation, the dis- 
tortion caused by coherent infall velocities takes a particularly simple form. On 
average, large scale fluctuations in redshift space 5s are enhanced with respect to 
real space 5 because of the radial velocity infall, so that 0/3 = (l+//3)5 

so that they are larger by a factor (1 +//3). This enhancement is anisotropic. In 
Fourier space: 

P,(k) = (l+/Ai|)2p(fc). (6) 

where P{k) is the power spectrum of density fluctuations 5, fi is the cosine of the 
angle between k and the line-of-sight, the subscript s indicates redshift space, and / 
is the velocity growth rate in linear theory. 

The correlation ^ [a, k) is related to the power spectrum by a Fourier transform: 

i{a,n) = JPs(k)e-^^^ (7) 

After integration in Eq.||7l, these linear distortions in Ps(k) produce a distinc- 
tively anisotropic ^{a,7z). At scales smaller than about 5QMpc/h there is a clear 
squashing in the correlation function caused by the peculiar velocity divergence 
9 field, this effect can be used to estimate /, for example by fitting the normal- 
ized quadrupole to the data (eg see ISl). For the sample with z=0. 15-0.30 we find 
/ = 0.48 — 0.83, which corresponds to Q,„ = 0.24 — 0.32 when we assume standard 
gravity (7= 0.55). 

Redshift distortions in the linear regime produce a lower amplitude and sharper 
baryon acoustic peak in the LOS than in the perpendicular direction because of the 
coherent infall into large scale overdensities. This is illustrated in the left panel of 
Fig|2] A characteristic feature of this effect is a valley of negative correlations (in 
blue) on scales between tt = 50 — 90 Mpc/h, which is in very good agreement with 
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Fig. 3 Signal-to-noise ratio 
in £,{a,Tt) for the z=0.15- 
0.30 slice. The color scheme 
denotes (S/N)^ multiplied 
by the sign of the signal i.e. 
negative values correspond to 
a negative signal. The triangle 
highlights the region 7t > a, 
which receives little weight 
in the monopole. The mean 
(S/N)^ in radial bins in this 
region is shown in Fig|4] 
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our measurements from real SDSS data. Such a valley is absent without redshift 
distortions. 

In Fig.|3] we show the signal-to-noise of ^ in the (7 — tt plane for the redshift slice 
z = 0.15 — 0.3. This complements the £, {a, n) signal plot in the right panel of Fig. [J] 
The signal-to-noise shown in Fig.[3]is for each pixel of size 5 Mpc/h by 5 Mpc/h (the 
same pixel size is used in Fig.|2]i. Note that there is covariance between pixels, and 
so this figure should be interpreted with some care (see |f3|). Nonetheless, it demon- 
strates the high quality detection of a BAO ring in the a — n plane. The triangle 
highlights the region n > G, which receives not much weight in the monopole, but 
where the BAO ring still shows up nicely. Note that the {S/N)^ shown is modulated 
by the sign of the signal: the (blue) valley of negative correlations at tt ~ 50 — 90 
Mpc/h - in accord with the predictions of the Kaiser effect - are detected with signif- 
icance as well. The overall coherent structure of a negative valley before a positive 
BAO peak (at just the right expected scales) is quite striking, and cannot be easily 
explained away by noise or systematic effects. 

The evidence for a BAO peak in the monopole is quite convincing (see |l2][3]|4l). 
The data follows the model prediction and produces a clear Qt, detection which oth- 
erwise (without the BAO peak) is degenerate with other cosmo logical parameters. 
But the monopole signal is dominated by pairs in the perpendicular direction a > n 
and here we would like to assess if the BAO peak is also significant in the radial 
direction. We do this by studying the signal-to-noise ratio in £, (cj, n) for n > a. In 
Figl3]this corresponds to the region inside the over-plotted triangle. We do the mean 
signal-to-noise inside the region tt > C7 as a function of the radius = + (7^, in 
radial shells d±ds of width ds = 2.5 Mpc/h: 

mean{S/N)^ ^ {Sign) {S /N[n,a]f (8) 

s±ds 

where {Sign) is the sign of the signal. When the signal is negative this gives a neg- 
ative contribution to the mean signal-to-noise square. If the signal is dominated by 
noise, positive and negative fluctuations will tend to cancel and reduced the mean 
{S/N)^. Results are shown in Fig|4] The mean signal-to-noise is always larger than 
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Fig. 4 Mean square signal-to- 
noise ratio averaged in radial 
bins with 71 > <J, shown as a 
function of the radial distance 
7t. Note that we include the 
sign of the signal in doing 
the mean, see Eq.d), which 
cancels the noise contribution 
(mean negative values are 
shown as a dotted line). 
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unity both in the negative valley between 50-90 Mpc/h and also around the BAO 
peak, where the mean (S/N)'^ approaches 2. This clearly indicates that the BAO 
peak is also significant in the radial direction and it also has the shape that is pre- 
dicted by the models, with a negative valley and a positive peak that extend in a 
coherent way over the expected lengths. It is unlikely that noise or systematic errors 
could reproduce these correlations. Similar results are found for the other redshift 
slices, with less significant detection for the middle slice. 
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